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ABSTRACT 

This report summarizes efforts to characterize the measurement of 
conductive mesh and smooth surfaces using proximity measurements for a 
dielectric resonator. The resonator operates in the HEM n mode and is 
shown to have an evanescent field behavior in the vicinity of the 
sample surface, raising some question to the validity of measurements 
requiring near normal incidence on the material. in addition, the slow 
radial field decay outside of the dielectric resonator validates the 
sensitivity to the planar supporting structure and potential radiation 
effects. Though these concerns become apparent along with the 

sensitivity to the gap between the dielectric and the material surface, 
the basic concept of the material measurement using dielectric 
resonators has been verified for useful comparison of material surface 
properties. The properties, particularly loss, may be obtained by 
monitoring the resonant frequency along with the resonator quality 
factor (Q), 3 dB bandwidth, or the midband transmission amplitude. 

Comparison must be made to known materials to extract the desired data. 
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CHAPTER 1 


INTRODUCTION 

1 . 1 INTRODUCTION 

This report summarizes the research effort to evaluate the potential for using a dielectric 
resonator in the vicinity of a planar, conducting, material sample to determine the loss properties 
of the sample. This introduction and resultant conclusions summarize the extensive work 
documented in the appendix and submitted as a thesis for the MS degree at Virginia Tech. The 
thesis develops the theory for approximate analysis of dielectric resonators based on a review of a 
variety of methods expounded in the literature. The technique of Itoh [4] provides useful result, 
but is found to be a bit too approximate for the problem at hand. A full problem development 
could have been pursued, but was not felt to be necessary for the justification goal of the report. 
In order to obtain useful results, a middle-of-the-road approach was taken which addressed some 
of the techniques suggested by Itoh and added a more complete description to the capped 
terminations of the resonator. In particular, the consideration of bringing a planar material 
within the proximity of a planar end of the resonator is included in the report. 

1.2 BASIC APPROACH 

To attack the problem at hand (see Fig. 2.1 of Appendix), the resonator was first 
considered to be a dielectric waveguide. Thus, the properties of a circular dielectric waveguide 
were developed for an infinite dielectric rod. The properties of the rod analysis were not fixed, 
but were developed for use in an iterative process as would be required for the determination of 
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the dielectric resonator properties, those of a truncated dielectric rod. Several classic cases were 
considered in order to provide validation to the code. The validation simply involved the 
determination of the fields and longitudinal propagation constant of the dielectric rod for different 
modes as expected from the literature. 

The next step in the development was the truncation of the dielectric rod to represent a 
circular dielectric resonator. The ends of the resonator were modeled as simple connections to 
open space. The boundary conditions at the ends were approximated by either the electric field 
or magnetic field mode dominance characteristic of hybrid dielectric waveguide modes. This 
boundary condition was thus represented by a reflection coefficient which was iteratively added to 
the dielectric rod waveguide delay to obtain a round-trip delay having a multiple of 360 phase 
shift. The iterative process consisted of searching for the complex frequency at which the proper 
phase shift occurred. In addition, the properties of the reflection coefficient magnitude and the 
dielectric rod propagation loss directly impact on the imaginary part of the complex frequency 
and thus the dielectric Q of the resonator. Thus the Q for a particular mode may be determined 
and used in an estimate of the loss. 

To obtain the material loss estimate, the reflection coefficient at the end of the resonator 
adjacent to the planar sample was modified to include the sample energy absorption. Thus the 
overall modification in the resonator Q was evaluated. Sample results for the analysis in the 
presence of a lossy material sample is included in Section 3.5 of the appendix. 
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CHAPTER 2 


CONCLUSIONS 


This report has summarized the potential for using a dielectric resonator in the vicinity of a 
planar, conducting, material sample to determine the loss properties of the sample. The theory 
for approximate analysis of dielectric resonators is reviewed, incorporating a variety of methods. 
The technique of Itoh [4] provides the most straight forward approach to the problem, but was 
found to be a bit too approximate for the problem at hand. Itoh’s approach was extended to 
address a more complete description with a truncated circular resonator. An additional extension 
included the ability to bring a planar material within the proximity of the planar end of the 
resonator. 

The resonator was first treated as an infinite dielectric waveguide. The properties of the 
rod analysis were iteratively used in a program along with an estimate of the reflection properties 
at the resonator ends to determine the basic complex frequency of the dielectric resonator, the 
imaginary part of the resonate frequency directly providing an estimate of the resonator Q, 
quality factor. The ends of the resonator were initially modeled as simple connections to open 
space. The boundary conditions at the ends were approximated by either the electric field or 
magnetic field mode dominance characteristic of hybrid dielectric waveguide modes. This 
boundary condition was thus represented by a reflection coefficient which was iteratively added to 
the dielectric rod waveguide delay to obtain a multiple of 360* phase shift. The iterative process 
consisted of searching for the complex frequency at which the proper phase shift occurred. The 
imaginary part of the complex frequency was used to determine the resonator Q and thus the 
resonator loss. 
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The properties of planar materials adjacent to the resonator end were used to evaluate 
the effects of the resonator loss estimate by modifying the reflection coefficient at the end of the 
resonator adjacent to the planar sample to include the sample energy absorption. Thus the 
overall modification in the resonator Q was evaluated. Sample results for the analysis in the 
presence of a lossy material sample indicate a change in the resonator frequency and Q, thus 
providing support for the use of the corresponding resonator transmission amplitude for 
comparison of material properties of materials set adjacent to the resonator end. It is important 
to recognize that the only information which is directly available is a comparison with other 
samples. Thus, several known materials must be used for comparison in order to calibrate the 
effective measurements used in the system. The data is found to be very sensitive to outside 
influences and sample placement distances. Extreme care must be taken in the placement and 
surrounding environment of the measurement system in order to obtain repeatable and useful 
results. If adequate precautions are taken, the results obtained from several material 

measurements will provide useful information on the quality of a mesh conductor used for the 
reflector of a space deployable antenna system. 
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DIELECTRIC RESONATOR IN THE PRESENCE 
OP A LOSSY CONDUCTOR 
by 

Scott B. Johnston 
Dr. William A. Davis, Chairman 
Electrical Engineering 
(ABSTRACT) 

This thesis develops a method for obtaining the complex resonant frequency of a post dielectric 
resonator in the presence of a lossy conductor. A full field analysis is performed on an infinite 
dielectric rod from which the complex propagation constant and modal solutions are found. 
Using a single dominant mode (HEM„), the boundary conditions at the end of the resonator are 
enforced, to obtain a complex reflection coefficient. Using the propagation constant from the 
infinite dielectric rod and the reflection coefficient derived from considering the dielectric-air 
interface at the resonator ends, a two dimensional search is performed to find the complex 
frequency for which the gain/phase criterion of the resonator is satisfied. In the final step, 
boundary conditions are enforced for a lossy conductor at a distance As from the dielectric which 
yields the objective - the complex resonate frequency of a post dielectric in close proximity of a 
lossy conductor. 
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Chapter 1 


Introduction 


It was recognized in the 1960’s that high-Q dielectrics held high promise in the area of 
microwave electronics. The advent of temperature-stable materials such as barium tetratitanate, 
allowed dielectrics to be used in microwave systems such as oscillators, filter banks, and 
multiplexers. Following the temperature-stability breakthrough numerous papers have been 
written on both material aspects of dielectrics and the electro-magnetic field models of dielectric 
resonators. Various scientists have devised schemes using high-Q post dielectric resonators 
between two conductive shields to measure physical properties such as dielectric constant and 
loss tangent of the dielectric at frequencies of interest [l;3]. The modeling work in this thesis 
will reverse the scheme by using a fully characterized dielectric to determine the electrical 
conductivity of a single conductor in close proximity of the post dielectric resonator. 

The work of this thesis is similar to work by other authors [4;-;7], in the attempt to 
derive the field distribution and to determine the resonant frequency of the post dielectric 
resonator. Itoh [4], considering the TE 0 j mode, makes several approximations about the field 
configuration that yield a simple numeric method for determining the fields in and around the 
post dielectric and the resulting resonant frequency. Due to the similarity of basic principles of 
Itoh’s work and this thesis, Chapter 2 reviews his paper, "New Method for Computing the 
Resonant Frequencies of Dielectric Resonators" [Itoh, 4]. In another paper, Marek Jaworski [7] 
achieves similar results to Itoh’s by applying a Green’s function in an 
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elegant, but rigorous mathematical model of the dielectric resonator. While Itoh’s and Jaworski's 
works address general characteristics of the dielectric resonator, this thesis considers a very 
specific problem: the complex resonant frequency of a post dielectric resonator of specific profile 
in close proximity of a lossy conductor. 

The goals of this thesis are threefold: first, derive the electro-magnetic fields associated 
with the specific experimental setup to facilitate a better understanding of how the evanescent 
fields couple to the surrounding experimental setup; second, develop a method of calculating the 
complex resonant frequency; and third, develop computer models to aid in the determination of 
the conductivity of the material being used as the ground plate. 

Chapter 2 presents background information pertinent to this work, introduces the 
resonator system under study, and outlines the method of development employed to achieve the 
three goals of this thesis. Chapter 3 develops the resonator system in tb~ four following stages: 
first, the field configuration for the infinite dielectric waveguide is derived; second, the fields and 
the resonant frequency for a free space resonator are determined; third, the fields and resonant 
frequency for the resonator in the presence of a perfect conductor are derived; and finally, the 
fields and resonant frequency for the resonator in the presence of a lossy conductor are derived. 
The final chapter provides conclusions and recommendations. 



Chapter 2 


Background 

2.1 Introduction of the Experiment 

The experimental setup under consideration is shown in Fig. 2.1a. The experiment has 
been devised to determine the electrical conductivity of the mesh used to shield the dielectric 
from above. The measurable parameters of the experiment are the S parameters of the system. 
Of interest in this experiment is the one port parameter s u which is return loss and the two 
port parameter s l2 which is insertion loss. At a resonant frequency, power will be coupled from 
one probe to the other. Resonance will be noted on the network analyzer when the return loss 
increases and the insertion loss decreases. The two probes of the network analyzer are 
positioned to couple lightly into the dielectric but not to couple significantly to each other [9j. 

The post dielectric rod is supported by rexolyte from below and the conducting material 
of interest in close proximity above. The rexolyte has a relative permittivity approximately 
equal to one and the thickness is assumed to be large compared to the penetration of the 
evanescent fields. The four conductor supports are far enough away from the dielectric not to 
couple into the radial evanescent fields. The validity of these assumptions will be verified in 
Section 3.3. 

Section 3.3 will show that the placement of the conductor should be within 0.5 cm of the 
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dielectric. One method of quantifying the mesh conductivity involves using several known 
conductor, as the ground plate. The mesh pressed directly against the dielectric provide, for 
easier analysis, but as will be seen in Section 3.4, is not necessary. What is necessary is a 
mechanism to provide a repeatable way to alternate conductors while maintaining the same air 
gap distance between the conductor and the dielectric. If the measurement for one conductor is 

made with zero air gap, the measurements for all other conductors should also be made with zero 
air gap 

As stated in Chapter 1, the first goal of this thesis is to derive the fields in and around 
the dielectric of Fig. 2.1a. The complexity of the derivation is greatly simplified by making the 
following assumptions: 

1) infinite conductor 

2) each region source free 

3) conductor and dielectric system in free space. 

These assumptions reduce the system under consideration to that of Figure 2.1b. Although the 
system of Fig. 2.1b is simpler than the experimental setup, the configuration of the system 
remains non- separable. In contrast to the infinite cylindrical metal waveguides, boundary 

conditions of the resonator must be satisfied at the cylindrical sides and at both ends. The 
problem of the non-separable geometry will be addressed in Chapter 3. 

Before reviewing the applicable literature, notation used in the remainder of the thesis is 
defined as follows: 

S - complex resonant frequency given by S = ( + jw 
v = resonant frequency 
< = damping factor 

k t = complex propagation constant given by = fi - ja 

a = attenuation constant 

* 0 = free space propagation constant 

Cj = complex relative permittivity of the dielectric 

a = conductivity of the conductor 
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2.2 Literature Review 

In the paper "New Method* for Computing the Resonant Frequencies of Dielectric 
Resonator" [4], Itoh makes several assumptions about the TE 0l electromagnetic fields of the post 
dielectric resonator. Itoh’s post dielectric resonator is shown in Fig. 2.2a. The assumptions 
made in Itoh’s paper are: 

1) majority of stored energy is in dielectric 

2) fields decay exponentially at the ends 

3) small amount of energy in fringe area, regions 5 and 6 of Fig. 2.2b. 

4) lossless dielectric (e t = e, - jO) 

5) lossless conductor («-*oo) at x = - I 

Itoh’s assumptions of a perfect dielectric and a perfect conductor result in a resonant frequency 
i s purely real. Using these assumptions, it is sufficient to match the fields that are tangent 
to the boundaries. Itoh assumes the fields in the four regions for an unknown s 0 are 


A j tinfi(z~z 0 ) J 0 (hp) 

region 1 

(2-la) 

A J sinfl(z-z 0 )K 0 (pp) 

region 2 

(2.1b) 

A 3 txp(-y(z-L)J 0 (hp) 

region 3 

(2.1c) 

Aitinhfcfz+iJJjhp) 

region 4 

(2 Id) 


where 

t* = »*-*? 
*• = w o^oK* 




Figure 2.2b Six regions of the resonator 
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The relative dielectric constant is and the thickness of the materiel below the resonator 
if , ^ fieldj in the region, five and rix are assumed to be »ero, or at least negligible. The 
remaining two components and H p are uniquely specified once H t is determined and are 

derived in terms of H t in Appendix A. Matching tangential component, at the dielectric 
boundary, performed in appendix B, leaves two eigenvalue equation. 

( 2 - 2 ») 

Whs) 

fiL = tan'flp) + coth(t) . ^ 2 - 2b ^ 

For the dielectric resonating in free space, t-*oo, 2.2b becomes 



(2.2c) 


Given a free space dielectric with height of L = 0.508 cm, radius of « = 0.63 cm, and 
^ = 40, numerical methods were used to solve the coupled Eq. 2.2a and 2.2c resulting in a 

resonant frequency of 


f t = 4.36 Ghz . 

Itoh’s work results in a quick method for determining the electromagnetic fields and 
resonant frequency of TE 0m modes. Itoh's method is not mathematically rigorous; however, his 
method for determining field, and resonant frequency is suitable for perturbation studies. 
Possible perturbation, include distance between ground plane and dielectric, change, in the 
dielectric constant, and change, in the physical profile of the dielectric. What the model does 
not provide i, a way to calculate hybrid mode, or to make provirion, for a conductor with finite 
conductivity. Method, addressing these two anomalies will be developed in Chapter 3. 
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2.3 PMC Dielectric Resonator 

In this section the free space resonator, shown in Fig. 2.3 a, is analyzed using the 
assumptions of a perfect magnetic conductor, PMC. This assumption has been a celebrated 
method used to determine electromagnetic fields, the propagation constant of dielectric 
waveguides, and the resonant frequency of dielectric resonators [10]. The PMC model assumes 
that the tangential magnetic fields are zero at the dielectric-air boundary. This simplifying 
assumption produces numerical solutions with errors as high as 25%. For example, it will be 
shown by a more rigorous method used in Section 2.2 and Chapter 3 that the PMC method gives 
a resonant frequency that is 20% lower than the actual TE on resonant frequency. The PMC 
method provides however a convenient "back of the envelope' 1 procedure for the quick 
determination of propagation constant and resonant frequency. 

The PMC derivation of TE^p modes begins with the solution of the Helmholtz's wave 
equation in a non-magnetic medium. The electric potentials that satisfy Helmholtz's equation in 
region one and region two (Fig 2.3b) are given by 

= A B l a (k^)(C t co,(k t z) * D x rittflkgt)) P* , o < « (2.3) 

am 

and 

*2 = eo.(k t z) + D z ,m(k t z)) , p > „ (2.4) 

am 

• * 

where B*' denotes a general Bessel function of the nth order. The subscript m denotes the mth 
solution for a particular n. The separation equation for the two regions are 

= (2-5) 

and 

*?=*£ + *?= m 0 «. 


( 2 . 6 ) 



Figure 2.3a PMC, Perfect Magnetic Conductor 


a 


< 

region § 1 

— >T 

region § 2 

dielectric 


free space 


Figure 2.3b Region one and two of the dielectric 
resonator 
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The TE field* are given by the following relationship* [10]: 


p - -J d 

(27a) 

H _ -i a * 

p jup dpdz 

(2.7d) 

jn d iO® 

^-Yp* 

(2.7b) 

H _ _-j_ a* * b 

n <P jupp d<pdz 

(2.7e) 

E,= 0 

(2.7c) 

H t = ^-A + 

x w dz* 

(2.7f) 


In region one, which is the dielectric, the fields must have finite values as p approaches 
zero. To fulfill this criterion the general Bessel function must be an ordinary Bessel function of 
the 1st kind of order n. The electric fields in the dielectric are found by operating on Eq. 2.3 
according to Eq. 2.7 which results in 


Ep - p P) (C i c0, (^i 1 ) + ,,n ( k z z )) ^ 

(2.8a) 

E\ = k x / n (k lfi p) (C l cot(k t t) + D x iin(k z z)) e>"* 

(2.8b) 

< 5 > 

II 

*5qT 

(2.8c) 

H\ = W < C 1 - D ' 

(2.8d) 

W> ( C t •'*(***) D 1 co, ( k z*» 

(2.8e) 

, l 

tfj = ^ J n (t lfi P) (Cj co,(k t z) + D, zinfk^)) 1”+, 

(2.8f) 


with the prime denoting the derivative with respect to the argument. 

In region two, which is free space, all fields must decay to zero as p tends towards 
infinity. Thus, in region two the Bessel function must be a modified Hankel function with an 
argument of j k^^p. The modified Hankel function with argument represent* an outward 

decaying wave provided that the argument is purely real. The field* in region two are found by 
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operating on Eq. 2.4 according to Eq. 2.7 which leave* 


fij = £v>V) (c z eo, M + D z «»(**)) j* V **) 

fi£ = ( c z co '( k * x) + D z 

fij = 0 (2.9c) 

Hj = KiCiV "<*W * D l W‘WV W < 3 -«> 

flj = a^L, K„(3i J< />, I-C, ,mfVi - 0* co.^wV"' <*•*•) 


and 


= ^ ^ coi^r; + D 2 ,in(k l z))<> n + 

The boundary condition as defined by the PMC model, at p = a, z = 0 and z=h are 


(2.9f) 


Hfa = •)=!%' = •) ~-0, 
H\(‘ =0)= = 0) = 0 , 

//Jr* =h)=H\(x=h) = 0. 


(2.10a) 

(2.10b) 


(2.10c) 


Enforcing boundary condition (2.10a) gives 

a^j W> «*<W + *>. «'*' = • 


which simplifies to 


J n(V> = 0 
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Define , = kj p , as the mth root of the Ordinary Bessel function Jn to give a solution 


1 



(2.11) 


Enforcing boundary conditions, (2.10b) and (2.10c) give 

c t = 0 

and 

D, ««n (kjh) - 0 

To obtain a non-trivial solution it is required that 

tm(k t h) = 0. 

Eq. 2.13 is satisfied if and only if 

k t h = px 

or 

p - 0,1, t, ... 


( 2 . 12 ) 


(2.13) 


(2.14) 


Substituting equations 2.11 and 2.14 into the separation equation (2.5) the resonant 
frequency is derived as 

*f = = "f Wc 


(uS, = sl<W + O 1 - 


(2.15) 


Where the indices n and p may be any non-negative integer and m may be any positive mteger. 
A s imil ar procedure is undertaken to find the resonant frequency of the TM^p modes (2.15). 



(2.16) 


„ \Thi , _ 1 _ 




<¥»' 


Th« symbol ^ is the argument for the sero’s of the derivatives of the Bessel function and the 
indices n and p may be non-negative and m any non-negative integer. 

The dielectric used in NASA’s experiment has a radius of 0.63 cm, a height of 0.51 cm, 
and a dielectric constant of 40. Using these dimensions, the resonant frequency and propagation 
constant have been determined for several modes. 


mode 

T^oio 
™„o 
TE ,n 
Win 
7*ou 


f t COuJ 
2.88 
4.59 
6 53 
5.14 
6.53 


U/m) 

o 

o 

615 

615 

615 


2.4 Complete Solution 


Before considering the details of the derivation and the approximation that allow for an 
accurate solution with finite time and resources, insight may be obtained by considering the full 
scale problem without approximation. The system consists of a lossy dielectric resonator having 
an excited hybrid mode operating in close proximity to a lossy conductor. Fig. 2.4a illustrates 
the difficulty. In contrast to an infinite waveguide structure, the dielectric resonator system has 
a non-separable geometry. The first step of the analysis is to break the system into ten regions 
as depicted in Figure 2.4b. The Helmholti’s equation is solved in each of the ten regions yielding 
ten sets of scaler potentials * aal Then, the fields must be matched at all interfaces. An infinite 
set of hybrid electro-magnetic fields are found for each region. 
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Figure 2.4a Experimental set up 
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Figure 2.4b Regions of the dielectric resonator 
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Region one of the dielectric resonator has an infinite number of mode., although only one or 
» few may dominate. A, in the other region., all but a few of the mode. repre.ent evane.cent 
mode,. Evanescent mode, that radiate into free .pace repre,ent stored energy and do not effect 
los. in the system. Evanescent modes that are terminated by a !o».y conductor may however 

have power transfer thus representing a loss mechanism. 

To complete the solution, mode-matching technique, are applied at the interface 
between ail regions. The boundary conditions state that the tangential component, must be 
continuou. at the region interface, except at the source interface. As an example, consider 

region 6. 


1) Fields must decay to zero as oo 

2) Field, must decay to zero as ?-*oo 

I 

3) Tangential components must be continuous with the tangential components 
of repon 6 and 7. 


The source* me, be creeled in phe*e end out of pbeee to determine the .render 
cherrclerietics of the re.on.tor ». en dl.rn.te procedure, the complex re-nen. frequent, of . 
*ource free repon me, be determined end the Q of th. re.on.tor computed b, the redo of the 
red to th. imefiner, frequent, term. The mode metchin, procedure outlined ebo.« doe. not 


guerentee re.ooence for the eource-free problem, rether th. procedure mud be repeeted until . 
frequent, i. found for which the field. ...id, dl of th. equdiem.. Thi. m*h. be dded in term, 
nf the bounce peth between the re.on.tor end, providin, . «ein, with the requirement thet the 
loop t em he of unit, moputud. nd . net phe,. then,, of r. Th. U»P «dn crit.n. for 
re.on.nce will be .xpntded in Ch.pter 3. It mud be noted thd the problem outlined .bone, 
dthoufh extremely ripnon., i. dill onl, d, drpr.ani.Uon to the red world problem. Two 
wiumpUon. which delin.de .hi. point d. the infinite pound plde did 4. point .ource. In 
order to work the problem in . finite dnoun. of time, dlditiond dmmpUon. d« necedd,. 
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Each assumption which facilitates easier mathematics must also be justified. In the following 
development the justification will be that the assumptions deviate slightly from the physical 
problem and introduce negligible error. 

2.5 Outline of the Development 

To facilitate a better understanding of the operation of the dielectric resonator system, a 
five-step "ground-up" development will be employed. The five steps of the development are the 
analysis of the 

1) infinite dielectric rod 

2) TE 0 , and TM 0 , resonator 

3) HEM,, mode resonator 

4) dielectric resonator in the presence of a perfect conductor and 

5) dielectric resonator in the presence of a lossy conductor. 

In step 1 a full field analysis of the infinite dielectric rod (IDR) will yield a method for 
determining the propagation constant i, for guided modes. As a result of the derivation, given a 
complex frequency, radius a, and a complex permittivity e d , the complex propagation constant 
may be obtained. Using results of step l, step 2 derives a formula for the complex reflection 
coefficient at the resonator’s ends. Using the reflection coefficient the criterion for resonance will 
be forced to yield a complex resonant frequency for the cases of TE 01 and TM 01 . Step 3 will use 
similar techniques of step 2 to determine the complex resonant frequency of a dielectric operating 
in free space at the lowest order mode hybrid HEM,,. Step 4 introduces a perfect conductor to 
the dielectric system operating with HEM,,. In the final step, analysis of a resonator system 
with a lossy conductor will yield a formula for an equivalent complex reflection coefficient. 
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Once the complex reflection coefficient if calculated, ftep* 2-5, the complex resonant 
frequency of the dielectric resonator can be computed. The method for computing the complex 
retonant frequency is based on loop gain and phase criteria. Similar to electric oscillator and 
optical laser cavities requirement# for resonance, the loop gain of the system must have a gain of 
1 and a change in phase of zero. Specifically for the dielectric resonator, loop gain and phase 
refer to the electro-magnetic fields, illustrated in Fig. 2.5a. The formula that enforces the 
criteria of gain and phase is given by 




e -w 


= 1.0 


(2-17) 


A two- dim ensional search will be used to determine the complex frequency, S, that 

■ninin.i t>. the error of Eq. 2.17. Each complex frequency yields a propagation constant. Given 
the propagation constant the error function is given by 

error = \ 1.0 - r„r l4 | . (2 18 > 


The two dimensional search is depicted in Fig. 2.5b. Once the complex frequency, S — ■< + ju 0 , 
is determined, the unloaded quality factor of the resonant system is given by 


< 

Programs have been written to simulate the operation of the dielectric resonator system 
with one of three modes exited under different physical conditions. Conditions of special interest 
are the perturbation of the dielectric loss mechanism the air gap distance between the 
dielectric and the conductor, and the conductivity of the ground plate. As a result of the work 
in Chapter 3, given the physical dimensions of the resonator system, the output of the simulation 
will be the complex propagation constant complex reflection coefficient, and complex resonant 


frequency 







Chapter 3 


Analytical Development 


3.1 Infinite Dielectric Rod 


TV (O.II of tbo fir.t it H 0 Of development are two: derive, without approximation, tho 
elec.ro.map.etic Bold, of tho infinite dioloctric rod and d«.i.« . method for determinin, tho 
complex prop**™. comtant t, Tho cylindrical dioloctric ...opud. uudo, c.u.id.r.«io« i. 
fbown in Fi|ure 3.1. Ropon 1 i. tho Infinite Dioloctric Rod (IDR) and repo. 2 i. tho mfimt, 
,p„, .urroundinp tho rod Tho d.ric.tioo bepn. b, min, ... .... of el.ctromap.<t,e potootiP. 
thu sali.fy Helmholu 1 . .caler equation in each of the two ropon.. Re.tnctint tho Bold, to finite 
tho .loclro-maptotic potential. mo mod to derive the hybrid Bold, for each of tho two 
ropon.. Applying the dioloclncdioloctnc boundary condition, th.t u»C«»tial Bold, mrnt b, 
continuou. .c... the bounds, remit, in . d* «> < «!«““»» and 5 unknown. Four of tho 
unknown, are Bold coefficient, A-D. »d tho fifth unknown i. the prop H .tion '»"*“* V The 
,„t two talk, of the ..Cion will b. to develop a method for rnl.in, 1, and ..l.in, the Bold 

coefficient. A, B and C relative to coefficient D. 

The derivation begin, with the magnetic and electric potential, that .atisfy Helmholtz-, 

equation in region 1 and 2. The.e potential, may be written in general form a. 


= AE^ x (k XpP ) e>"* 

am 

nm 


(3.1a) 

(3.1b) 
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region H 2 


region §] 



y 


Figure 3.1 Infinite Dielectric Rod 


am 
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(3.1c) 


(3.1d) 


The oot.lion fl» denote. . pnerel (unction of the otb order. The ntb.cript m denote, 

the n*h root of », .oln.ion (or . puticulu n. To be explicit tbe pr.pH.tion c.n.t»t 
repre. online mode .. in repon 1 «ould be .ritten l„„ To >.ti.fy bounds condition. .« .be 
dielectric -dielectric interface for HI it i. required th.t Thu., the propH.ti.n 

constant will be shortened to ^ 

A direct result of solving the wave equation is the separation equation. The separation 
equations for the two regions are given by 


,t _ t 2 f _ t * 

*1 - * 0*1 - * 1 ? 


+ tf — S/'oVo 


and 




( 3 . 3 .) 


(3.2b) 


For cylindricl ....pride., mode. tb.t exhibit .npriu .ufation — be pnr. TE or TM 
■node.. The mode, th.t exhibit u.priu .ufation the cylindricl dielectric ». (HEM) hybrid, 
commonly referred to u either HE or EH mode. HE i. o..d «h.n TE (H) mode, ore 
predominont mid EH (E) .hen TM mode, me predominmit. For dielectric ....pride, the 
domino mode ,. HEM,, foiled b, TE„„ HEM,,, u>d TM,, (.»]. TE mid TM mode, me 
.pecio] cm.., n=0, of the infinite hybrid .«t. Thu., con.ider.lion of tbe hybrid .«t fake. cue of 
ril po.aible mode. The hybrid mode, cm be .ritten fa . .uperpo.it, on of properly oper.t.d on 
r cod The hybrid electric ond mnnetic field, u. pven by 




H r 


- 1 3 #M 
~ p H 


+ 




(3.3b) 


(3.3a) 
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= + k*‘ (33c) "* = ' + w 4s < 33d > 

E, = {3M) ' ?) * t < 33( > 

In region I, the fields must have a finite value as p approaches zero. Thus, the general 
Bessel function is an ordinary Bessel function of the first kind order n. The electric fields in 
region 1 are thus given as 



(34*) 


(3.4b) 

*! = x ^ VV; S*< A ‘ 

(3.4c) 


where if, - if - if 

The same method is applied to yield the magnetic field components of region one. Eq. 
3.3 operating on 3.1 yield 


H l p =lA& SrnttuS) 1 <’"*<*'* 

(3.4d) 

»i = M ‘i» /.<*,✓,» * *’■*«*' 

(3.4c) 


(3.4f) 


The dielectric of region two is free space with e 0 and ji 0 and assumed to be infinite in 
the p direction. All fields must be finite as p tends toward infinity. This requires the general 
Bessel function in region 2 to be a modified Hankel function with an argument of The 
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modified Hankel function 
expansion* are given a* 


„,lh T«um-nt repreeentl » d.c.,to, «*»«• The field 


£> = / c K. (tos> * D '$ K ' (i ^ > 1 **' ' l, ‘ 

(3.5a) 

g} = (Clg} K,(,t- Ir f) > D(jt»K 0*tS)l 

(3.5b) 


(3.5c) 


wbere ^ _ L l . t l A* before, the magnetic field components in repoo two become 


wj = /c^ k,(jV^ fD * 0 f*. *> (i ^ } 1 * e 

(3.5d) 

h\ = ic (■-*,) tf.Ohs) * D ^e 1 

(3.5e) 

4 *o%k **,)'*'*•■ 

(3.50 

The held, den.ed obo.e represent on infinite number of pnmib!. field .nlu.inn, to 

addition to the .-dependence, the tangential components must be continuous 
dielectric boundary. These boundary conditions may be wntten as 

at the dielectric- 

1) H\ 6>=«, 0 < * < »,*)= H l (^ a ‘ 0 - * - *'■ 

(3.6a) 

£ t ‘ ©<#<*», O = ^ fr =a - 0 * ♦ - *’ ' ) 

(3.6b) 

y, < * < *.0 = H\(P=*’ 0 ^ ♦ “ ^ 

(3.6c) 

and 

*; E\(p=*,0< * < = %(?=*> 0 * * - 

(36d) 
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Substituting the field expansions into the boundary conditions, we obtain the following four 
equations: 

B if, = D if, (3 7a) 

» 4 ^ 4 (3 7b) 

A k ip^ »( k \p d ) + *h ~ C jhprf mtitip*) + K*Ghp a ) (3.7c) 

&nd 

* iirejk J J k \» a ) f B k \p* J k \p a ) ~ c KGhp*) + D fop* ' u 0*2f*l (3 7d) 

These four equations can be set equal to zero and written in matrix form 

*_«=» < 38 > 

The quantity x is the eigenvector of the coefficients given by 

A 

B 

x- (39) 

C 

D 


and matrix Z .„ is given by 
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0 


0 

ApW^p*) 


0 


0 

k if>fy k \p a ) 


■Jh pKtihp*) 



k \A( k \„ a ) 

0 0 

■it 


where the matrix elements are a function of i r The k p component for regions one and two are 
given by 


k 






&nd 




= N " 2 


(3.10a) 


(3.10b) 


The subscripts nm are added to emphasize that there are multiple solutions 
corresponding to the nth harmonic in <f> and mth root of the matrix equation. There are two 
possible solution forms for eigensystem (3.8). The first trivial possibility is that the eigenvector 
is equal to zero. The second possibility is a solution requiring the determinant of Z equal to zero. 
Since 3.8 is a function of t,, to have a non-trivial solution there must exist a particular A, that 
forces the determinant of Z to zero. The technique used to find the particular complex 
propagation constant A, is outlined below. 

A direct search algorithm is used to find the k t root by first searching the dominant 
phase term for a minimum determinant and setting 


1| twin • 
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The method used to find 0 searches starting with /},. increment* by A/S, and end* at /Sphere 


and 

^1 = 4 W f# 1 O e o e T 


Any 0 outside this specified range will be imaginary, representing an evanescent wave. As 0 is 
incremented, the magnitude of the determinant Z Bm , denoted by | Z,^ | , is calculated to 
determine which fi yields the minimum determinant magnitude. The results of this method is 
shown in Pig. 3.2 for 

freq = S. 54 Ghz 
a — 0.6S5 cm 

and 

e t = 40 -)0 04. 

Once dmin established, A, is set equal to l t = 0 m ,„ and a second search is 

undertaken to find the real part of the complex propagation constant. The real part of A, is 
incremented in search of the particular A, that minimizes | Z, m | (Fig. 3.2b). Once the two 
dimensional search is completed the propagation constant is given by 


(3.11) 
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The retulti of the two dimensional search for this example provides the solution 
k t - 319.0 - jO. 336 for the HEM n mode. 

The negative imaginary component implies attenuation of the forward propagating 
waves. Having computed the propagation constant k 2 , the p component can be determined using 
Equation 3.10. The p components are 

= 344.0 - jO. 00837 
k^p = -0.345 - j0.336. 

Once is computed, it is insightful to determine the coefficients A, B, and C relative to 
D. The coefficients are related to the magnitude of the field components and can be determined 
to within a constant without having information about the source. With knowledge of the 
source, a Green's function technique may be used to determine the four coefficients. For the 
purpose of this paper, it is sufficient to calculate the coefficients A, B, and C relative to 
coefficient D. 

To determine the relative coefficients, matrix system (3.8) is written in the form 


0 

f\7 

0 

At 

h\ 

0 

/j3 

0 

h\ 

hi 

As 

A 4 

f<\ 

fa 

A 3 

At 


(3.12) 


which hM the form 
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Zx =0 


Taking 


the LU decomposition of Z leave, the matrix in the form 


- 


- — 



?n / u 


>4 


0 

0 lii h 3 ^34 


fi 


0 

0 0 /j3 /34 


c 


0 

0 0 0 fa 


D 


0 

- 


- 


L -1 


(3.13) 


(3.14) 


The last equation of the matrix system is 
D*f u = 0 , 

where Af „ m 0. Eq. 3.14 allow, D to be assigned any value. The simple choice of D = 1 may be 
used to obtain a 3x3 matrix in the form of 


ll =*, 


(3.15) 


-1 


- 


— 

?H / J 3 /l 3 


A 


-f 14 

0 ?J 3 A 33 


B 

= 

" / 3 4 

0 0 ?33 


C 


"/ 34 



- 


__ 


(3.16) 


where , U . rente*, end .he ,..«r * me. Urn* » >W A. B. C reletive to D. Cedin, « A. 
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work outlined above Allows the field components of equations 3.4 and 3.6 to be plotted versus p. 
Usinf the dimensions and properties of the NASA experiment, the 6 field components versus p 
are plotted in Fig. 3.3 for the HEM„ mode. So that the electric and magnetic field strengths 
can be compared, the fields are given in units of Several observations can be made to 
intuitively check the results including 

1) tightly bound fields 

2) continuity of tangential fields 

3) discontinuity of the E p components by a factor of e r . 

4) predominately TE (no E J 

The field pattern for HEM n is shown in Figure 3 4. The knowledge of the field configuration for 
the specific dimension of the experiment is essential in the further development of the dielectric 


resonator system. 
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3.2 Resonator Method for TEjj 


A> lttttd „ sect,.. », d. *- “ *— ^ ph,,,Cd dta ™ i0 “ 

„ d pmpetie. e. th. dielectric. * — . - ,UCh “ *• “ 

th. pr.be., P u,. - — c. - — -* — " ° D ' ** 

. . t u, TE mode. In thi* aection, the complex 

. .l A h» dominant i* tne m>oi IUOU 

low order modea that may be aomw«u 

„ ffl ,,ti.» — * b. derived for . free dielectric (Fi,. 3.5) M 

damirur* TB„ field confiierelioa - in Section 1 . 1 . 

B, eeneiderioi both forwerd end reflected ..... - *• — — -*• * = * 
exprereion fer Ur. refiecrie. coefficient. F... be forint. *• — “ S "“°° ^ ** 
reflection efficient -iU be n«d .. Uererivei, for *. *— * «— « 

frequency. The derive*., be*, b, u.in, the conic from the infinite dielectric rod of Set... 

, . The teeth order electro-m^etic field. ere — *—*-••* * “* 

C „ m, ecpiel ter. .hieh re.nlt. in pne — electric ...... ^ TE., field. In re*.. 

one are riven by 


and in region two by 

h 1 , = 


(3.17a) 

(3.17b) 

(3.17c) 

(3.18a) 

(3.18b) 

(3.18c) 
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Alon( the dielectric-air interface, p= a, the tangential components from region one to region two 
must be continuous. These two boundary conditions are stated in by Eq. 3.19a and 3.19b. 

I) E l /p=a) = B^(p=a) 
t) H\(p=a) = H\(p=a) 

Dividing (3.19b) by (3.19a) results in 

Vogt g*) =0 

or 

oK*) + hpKtdhp*)* \K«) = 0 (3.20) 

It should be noted that Eq. 3.20 is a special case of the eigensystem (3.8). If the matrix 
approach is used, setting n-0 would result in two possible solutions corresponding to TE^ and 
TM 0m . Satisfying Eq. 3.20 is an infinite set of eigenvalues ^^with the subscript m taking on 
values of 1 to infinity. The corresponding k components can be found from 

”Om 

(3.21a) 
(3.21b) 

From this point , only the first solution, m=l, will be considered and the subscript m=l will 
be implicit. The development proceeds by matching field components at the end of the 
resonator shown in Fig. 3.6. At interface, 1 =^, the tangential component* of region one mutt 
equal the tangential component* of region four. The fields in region one are composed of a 
forward and reflected traveling waves. The reflected waves are given by substituting i, - -i, 
into equation 3.17. In region four it is assumed that the fields have the tame component as 



(3.19a) 

(3.19b) 
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U.. field* in repon «... Thi. -eumptio. in *«*ct *««. th« In.,* field*, depicted in Fi«. 3.«. 
o, the reeonetor. The errer u^.etcd «i«h thi. ».«n.p.i.n .dl be «ndl il the Irin,. ere. i. 
ln4l | comp „cd to the totel luriece ere. ef the monetor end end the eeen.ee.nt field, dee., 
repidl, in the . direction It foil... thet = repon (.or, the propH-ion eon«.nt 
is obtained from 


n iT 

*4» — - \ w o^o*o l fi 


(3.22) 


The ne.eti.e imopnor, toot i. eho«. to pee dee., *»■«»*•“ The bounder, 

condition for the TE 01 mode at the end interface ia given by 


J )£?(*=$) + = 

v **?(*=%) + H U i= i) = fy i= h) 


(3.23a) 


(3.23b) 


The corresponding reflection coefficient 


I'm 


is defined at i — 


L 

2 


aa 


r,« = 



(3.24) 


Uaing i = 0 aa a phaae reference Eq. 3.23 can be rewritten uaing Eq. 


3.24 which yielda 


fiJ + r^r 14 ; = E} 

i 


(3.25a) 


(3.25b) 


Using *=Jf aa the phaae reference, the field components are 
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K* = »i Mis> 

e} = ’ 

’ ilM> 


and 


u* — a j> n B 

"J “ fl « w M 0 'o <*/// 4 


The boundary conditions (3.25a) can be rewritten as 


Bkj/ofk^Xl + r u J = B t k lfi / 0 (k lpP ) 
B(i + r u ) = B i . 


Boundary condition 3.25b can be rewritten aa 


b^t 4<V> - r„; = W; 




Dividing 3.27 by 3.28 yields 

£^£jli = _i_ 

V* - r l4 ; * 4I 


r 

,4 ~£+C 


^o,j 


(3.26a) 

(3.26b) 

(3.26c) 


(3.26d) 


(3.27) 


(3.28) 


(3.29) 


The reflection coefficient for a single dominant TMqj can be derived using Eq. 3.4 and 
setting n, B and D equal to zero. Proceeding in similar manner , as for the reflection coefficient 
for transverse electric, the transverse magnetic reflection coefficient is 


r !4 


= e £^ll ^ 

«r*, + *. 


(TM 0X ) 


(3.30) 
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H..HX derived tb. complex refaction coefficient olio., the ripirithm developed id Section 
J.« to be implemented end jield. tbe complex re.onont frequencj. Shooin beloer ie the remit of 
tb. ,ei»/phue criterion W". «itb Itob’e method [t|. Tbe r— b« . dielectric content 
= t« - jO.Ot (lorn untent 0.001), 0 rodiu. of 0 «3 cm ond o beifbt L of 0 51 cm. 


mode 

method 

*1 

TE 0 »i 

PMC 

615.7 

TE 01 , 

Itoh 

354 .672 - jO. 353 

TE 0U 

Pha*e/Gain 

354.670 - jO.351 


S (g+W») Hi 
-2,737,800 + j2* 5.476xl0 9 
-2,200,000 +j2* 4.375xl0 9 
-2,200,000 + j2* 4.375xl0 9 


A. .hown in thi. .ection, computin* the complex re.on«t frequency for the dielectric 
runout mode TE„ cm be don. urin, Grin/Pbue enteric. CcSdcc md vriidit, of the 
Gein/Pheoe method i. pveo b, concurrence .itb lub'e method, .bicb i. ‘bn.reti.ril, 
equivrient. to th. next eection, otunlion i. turned to tb. prritlem rf computin, tb, reoonmri 

frequency for hybrid model. 


3.3 Hybrid Resonator 

A. mentioned in faction 3.1 tb. lonet order mod, lor . clindricri dteketri. ••»«*■*. 
i, ri» bjbrid HEM,, mode. Other nomenclriur, u..d to demrib, the bd order mod, i. HE,, 
«bicb denote, tbot th, mode i. dominontl, trevem electric to thi. mctic, tb. problem, 
eraocieted .itb aolvinf for tb. complex reflection coefficient lor hjbrid mode. «e ctlbud. 
Severri .ppr.oche. oddretin, the* problem. or, pre«n«d. Tb. problem for thi. .peci*. 
remnrior i. remixed b, eppl,in( tb. i»*ht «rin,d in tb, d.vekpmc. ri the inftoH. dielectric 


waveguide in Section 3.1 to the dielectric re*onator. 
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The task of formulating a closed-form solution for the reflection coefficient, as derived 
for the TE and TM cases, will prove to be formidable. In. fact, it will be shown that using a 
mode-matching technique, as employed for the TE and TM cases, will lead to a numeric 
ambiguity of four equations and three unknowns. Several approaches are considered to resolve 
this problem. The final approach will draw from information obtained in the study of the 
infinite dielectric waveguide in Section 3.1, and the field distribution unique to the physical 
dimensions and material parameters of the experiment under study. 

The derivation begins by matching the forward and reflected waves of region one (Fig. 
3.7) to the fields of region four. The boundary condition equalities may be written as 


E'S + B l p - = 

K 

(3 31a) 

EY + £}* = 

K 

(3.31b) 

h) + + H l ; = 


(3.31c) 

hY + Hi' = 

< 

(3.31d) 


As in Section 3.2, the fringe fields are ignored and a tightly bounded field configuration is 
assumed. The propagation constant in region four is 






(3.32) 


Using the electromagnetic fields of Eq. 3.4 and noting that the reflected waves have propagation 
constant -k t , the boundary conditions may be rewritten as 
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-^(*1 />?) + (> P))( e ** 3 + ^14 e " 


m 4 ^ 7 t ^ lpP )+ B i 4 y ^ v ^*'* 4 * 3 


f + **., W ^ € * 3 + ^14 e * 3 ) 


M^Sji / a (*,W + fl 4*ip4(*i^ ‘'* 41 * 

M ? 4 (* 1^^' ;4 ‘ 3 - r >4 M = 

M 4 ?^(* iW + * 4 ^ . 

^V^*' 3 - r M e * 3 > = 

(-A 4 ii^^,(* lp />^ + ® 4 W 0 /igp ^ 


(3.33a) 


(3.33b) 


(3.33c) 


(3.33d) 


Using the numeric methods developed in Section 3.1, the coefficients A, B, C, and D can 
be found relative to D for the forward and reflected waves. Knowing the propagation constant 
i 4| and the four coefficients, A, B, C and D leaves the four equations with the three unknowns 
r l4 , A 4 , B 4 Rearranging equation 3.33 such that the unknowns are on the left of the equality 
leaves 


( a %£ tiw * «? w^* s+ && <«v>^ 




+ B, Cf Wx? )>’'"*) - /»<V- , + ®£ W\S»‘ nt 


M 


(3.34a) 


~ r u 
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\i (*% r ^ }) ^ + 

+ ; n(V>» + ^ 


and 


r !4^ *1^ + } tS k \pP)) ^> f ^4^l/.(*l^> J U ^ 

+ ®4 J (i^p ^ n (^ipPJJ e ^ = + ®w 0 ^i Ai(*lp^^ e 


Bach of the constants is assigned a function designator to indicate the position such as 

A l (p) ~ ( A ij^ 4(*lW + • 


This allows for 3.34 to be written in the matrix form 
f\\(p) f\j(P) A 3 (p) 

f'ii(p) f'i'i(p) A z(p) 

Ai (p) fn(p) h*(p) 

U\(p) A i(p) U$(p) 





§\(p) 







9i(p) 


^3 


h(p) 


*3 


u (p) 


- - 





which is in the form of 


(3.34b) 


(3.34c) 


(3.34d) 


(3.35) 


Up) ■ = <J<W. 


(3.36) 


Several problems must be addressed. First, the problem is over specified by having one 
more equation than unknown. Second, the 4x3 system is a function of p. A decision must be 
made where at the resonator ends to match the fields. The restrictions of the choice are 0 < p 
' *, where p is the variable that must be chosen and a is the radius of the resonator. One 
approach is to rank the four boundary equations (3.31) in order of their importance and the least 
important equation could be "thrown out". The information obtained from the infinite dielectric 
waveguide could be used to make a judicious decision. A second approach using all four 
equations employs a more elegant method of a weighted integration over the radius. 

In this method the first step is to forward multiply the system of (3.36) by its own 
transpose. This step may be written as 

& (p) ' Z(p) •= &(p) G(p) (3.37) 

As a result of the multiplication the 4x3 matrix Z(p) of (3.36) becomes a 3x3 system (3.37). 
Next, to resolve the problem of which point p to pick along the radius, the equation may be 
integrated along p leaving 


| J*(p)&(p)' Z(p) “/ = \w(p)Sf(p) G(p) dfi (3.38) 

The function u/fr) is a weighting factor that is chosen to facilitate the integration. A judicious 
choice of w(p) could yield a reasonable closed form solution of (3.36). More realistically, numeric 
integration techniques would be applied to solve for the vector a. Integrating and realizing a 
solution is beyond the scope of this work. 

An alternate approach uses the insight gained from the earlier study of the infinite 
dielectric waveguide of Section 3.1. Pig. 3.3 shows that the transverse electric field components 
of the HEM|| dominate the components in the direction of propagation. Thus, in this approach 
the calculation of the complex reflection coefficient will be done by assuming HEM^ m boundary 
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condition. along the resonator, p = a, *nd TE,^ bounder, condition, at the resonators end.. 
The reflection coefficient for the TE ewe derived in Section 3.2 i. 

r - k * ~ ill (3-39) 

r H ~ k t + k u 

Implementing the algorithm developed in Chapter 2 u.ing Eq. 3.39 yield. Figure 3.8. The graph, 
.how the first two resonant frequencies of HEM n . The first two resonant frequencies are 
designated ff m and f ll3 The properties and dimensions of the dielectric resonator are 
e, = 40 - j0.04 
a -■ 0.63 cm 

and 

L = 0.51 cm, 

yielding the complex resonant frequency 

S = - 1.2 10 7 + j2* 3.535 10 9 (HEM U ), 

field parameters 

/_ = 3.54 Ghi 

k^ = 74- 1 J + JO 04 m' 1 

k Xp 344.15 i-jO.044 m' 

= -0.0S-jS09.86 m' 1 
i, = S18.8-j0.08 ro* 1 
l 3i - .0.036 +J3S6.07 m' 1 

and 

r 14 = -5.48- 10'* - j0.9977. 


The second order resonant frequency i* 
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S - -1.2 l0 7 +j2* « 44 10 9 (HEM,,). 

From Fig. 3.9a-d it is evident that a percentage of the fields are external to the 
dielectric. These evanescent fields in free space represent stored energy. From examination of 
Fig. 3.9, it can be determined that the external fields approach aero as p approaches 2 cm. Fie- 
3.9d displays how the fields decay at the resonators ends for a resonator in free space, suggesting 
the importance of the external fields is of paramount importance to the experiment under study. 
External fields must be understood in order to couple energy effectively into the system from the 
probes, not perturb the experiment with the experimental support structures, and perturb the 
fields in a repeatable way with the different conductors. Section 3.4 exams the effect of bringing 
a perfect conductor within close proximity to the resonators end. From Fig. 3.9d it can be seen 
that close proximity is within 0.5 cm. 
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3.4 Dielectric in the Presence of a Perfect Conductor 

ThU section focuses on the dielectric resonator in the presence of a perfect pound plate 
a. depicted in Figure 3.10a. It is the goal of this section to derive a method of computing the 
complex resonant frequency a. a function of the distance between the dielectric and the perfect 
conductor. The distance between dielectric and conductor will be referred to as As. Before 
deriving an expression for the reflection coefficient, it is useful to anticipate the results that the 

derivation will yield 

1) as As— ►oo r u = r i3 (r,j - free space resonator) 

and 

2) as |r l4 | = 1 (perfect conductor pressed against dielectric). 

The first statement indicates that when the perfect ground-plate is placed infinitely far from the 
dielectric, the reflection coefficient r u , becomes that of free space and yields the same resonant 
frequencies as those in Section 3.3. The second statement describes perfect reflection as the 
lossless pound-plate is pressed against the dielectric resonator. In the second case the dielectric 
pressed against the conductor should yield the same resonant frequency with the same transverse 
field distribution as a dielectric in free space with a length 2L (Fig. 310b). That resonant 
frequency for the free space resonant dielectric with length 2L will be the hybrid mode HEM na . 

The derivation for the reflection coefficient begins by considering the electric potential in 

region 4 

#;= + fl * M (3 40) 

One. eg**n the fringe field, ere ignore end it i. Meomed thet the '• ere equel in region, on. 
and four. Thui i 4g can be given by 
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Figure 3.10o Dielectric resonotor in dose proximity 
to a perfect dielectric 


Perfect 

Conductor 




Fiqure 3.10b Configurations yield equivalent resonate 
frequency and transverse field distribution 
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*4*- + W* r * ^ ( 341 ) 

Th« positive root is chosen to give decay in the -i direction. The ^ component of the electric 
field is given by 

£} = ^ )) + B * ' +A * »} ■ 

At the air-conductor interface, *= - ( jj+A* ), 

= ^/o( k \fiP)( A i co, ( k u(°)) + B i tin ( k u(°))J = 0 

or 

>»4 = * 

Using this result, the potential and Eq. 3.3, the TE 0 , fields for region four are given by 


fij = B t »*&«(*+%+** )) ( 3 - 42 *) 

K = ». -w w; "•^..<''4 w < 3 « b > 

«! = 8, Wlff) )) ■ l J ne > 

Using i = as a phase reference, the H p ’» and ’» in region one at i = ^ are 

Ej* = B ^/ofrxfS) ( 3 43 ‘) E V ~ r H fl (3 43b) 

H\‘ = B ^ ( 3 ,43c ) = - r i4 B sg£ • ( 3 43d ) 


where the reflection coefficient is defined *s 
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r - E i 


At the air-dielectric interface, s = the boundary conditions are 


ey + ej- = 


and 


h\+ + H l ; = h* 

p p p 


Substituting (3.43) , (3 44) and (3.42) and into the boundary equations yields 


B (T 14 + 1) = fl 4 iin(i 4l As; 


and 


Bk t (' r u + 1 ) = -i B i k u cot ( k u* 1 ) 


(3 44) 


(3.45a) 

(3.45b) 


(3.46) 
(3 47) 


Dividing (3.46) by (3.47) and solving for the reflection coefficient gives 


p = *i <aw (*4i A ^ (3.48) 

14 t, +>t 4l 


As As approaches zero, the reflection coefficient becomes 


hm T 

A z->0 


U ~ 






Providing that Eq. 3.41 yields an imaginary t 4l » as the perfect conductor is moved an infinite 
distance from the dielectric the reflection coefficient is 


hm 

As-*oo 


r l4 


_ k l(-j) j k it _ * ^ 4 » 

k x('f) - i k u k i~*~^ti 


As Az— »oo the reflection coefficient becomes identical to that of the free space resonator in 
section 3.2. Having solved for T 14 in terms of Az allows the plotting of complex resonant 
frequency of the first two HEM|j modes vs. Az as shown in Fig. 3.11 and Fig. 4.12 respectively. 

In the final stage of development Section 3.5 considers the dielectric resonator in the 


presence of a lossy conductor. 
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Figure 3.11 Perfect conductor at As from dielectric for f m 
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Figure 3.12 Perfect conductor at Az from dielectric for 
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3.5 Dielectric in the Presence of a Lossy Conductor 

Th. development in thi. leet .ection i. the ..-nee of the them: the poet dielectric 

ope-in. i« the pre.«»c of • >»"* *“ 

ln ,„ , ,h. rn.themet.cel tool, thet mop the then,, of the "»■» Muenc, due to th, 

chen|e 1. th. co.ductl.it, .1 the conductor The objective require, th, dorivotion of en 
equivelent reflection coefficient .hot inc.rp.r-e. the eir «d the h.», eductor (Fi,. 
3.13.1. To derive on expre«i.n for T.,, tren.mi».io. line them, 1. emplojed. Th. equivolent 
tren.mir.ion line problem 1. .bo.n in Fi«. 3 13b. Th. volt.,, mure. » trnminoUd into . 
tren.mii.ion line of cher.aen.ttc impede... end l..«tb e. Srttum «K of the t™-» 
line terminete. into . line »i‘h cherecten.uc impede... e, « I.nfth », uhich in turn b 

terminated into load Zf 

To formulete . .olution for the equivelent reflection coefficient, r„, multiple reflection. 
mu .t be eccounted for. Th. too reflection coefficient, ere 


r, 4 - reflection from the dielectric-air boundary 

I 

r 4k - r*fl* ction form th< 4ir ' ,0 ® sy conductor bound&ry 


Firet, it .hould be noted. T„ end To, ere function, of e rin|l« inUrf-e. For exmnple, 
F|4 , developed in Section 3.4, 1. unique end -pecifie to th, dielectric-rir mterfec. Second, T* 
Jli^e end .pecific to th. eir-conductor interfere. To further compile-. th, problem there 
„ infinite number of r.n.ctio». in the eir «ep. or in S««« 3 of the e»el.,.u. trmrini.rio. 
Th. equivelent reflect, on coefficient murt tek. in to -count ,«h of the reflection. 


mentioned above. 

Field theor, c .ol», the problem b, u.i», en mfinit, -Tie. to »lv. for the reflection, 
of th. eir ,ep. A more convenient epproech emplo,. • t-hmqu, urtd in control theor,. 



Figure 3.13b Transmission line and equivalent 
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Treating the air gap as a feedback mechanism, the equivalent reflection coefficient may be 
written aa 


eq 


_ r S-li - — 

— I 14 * .4^ Aa 


I + e 




(3.49) 


Which reduces 3.50 to 


r +r /^4 « Ai 
r - E iij ± ik e 

l eq~ 


s i r p 
* + 1 14 1 4k e 


( 3 . 50 ) 


The reflection coefficient r i4 , developed in Section 3.3 is 


r»« = 




-k 


ii 


■*" ^4l 


( 3 . 51 ) 


The reflection coefficient for the air-conductor boundary is 


where k^ s is 

k 4l = - - k ] C 

Care must be taken in the choosing of the proper root. The negative sign is chosen to give decay 

to fields propagating in the positive z direction. 

As in Section 3.2, the fringe fields are ignored and tightly bounded fields are assumed. It 
will also be assumed that the rho components of the electric field in the air, lossy conductor, and 
dielectric are equal. The lossy conductor is assumed to be a good conductor having the 
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properties [14] 


^ > tOOO 

and 

—f~ +} T~T~ 



(3.53) 


Using the newly developed propagation constant and equivalent reflection coefficient, the 
complex resonant frequency can be calculated for the conductor at different conductivities. Pig. 
3.14 and 3.15 are the result of using the following parameters 


^ = 40 - j0.04 
As = 0.2 cm 
a = 0.635 cm 


and 


A = 0.51 cm. 









Chapter 4 

Summary and Conclusions 


By deriving the electro-magnetic fields in and around the dielectric and establishing the 
linkage between the complex resonant frequency and the parameters associated with the 
dielectric resonator, a thorough understanding of NASA's experiment has been achieved. The 
derivation and illustration of the manner in which the fields propagate in the specific dielectric 
under study; the strength of the evanescent fields; the effect of the conductor on the resonant 
frequency; and how conductivity of the conductor effects the Q of the system, have given useful 
insight to the experiment. 

The electro-magnetic fields were determined by enforcing boundary conditions of the 
resonator system. Section 3 3 showed that the fields at the resonators end's are predominately 
TE for the dielectric under consideration. This characteristic allowed for the reflection 
coefficient for the TE mode to be incorporated into the determination of the complex resonant 
frequency of the hybrid mode HEM,,,. This determination of the relative field strength of each 
component has provided an understanding of how the fields propagate and are distributed m and 
around the dielectric. Section 3.3 also showed that at the dielectric air interface the angle of 
incidence was approximately 75*, as measured from the normal. 

Through this analysis, two important details regarding the design and application of the 
NASA experiment were revealed from the knowledge gained by deriving the evanescent fields. 
First, the strength of the radial evanescent field, (pg 50) determines the distance the support 
structure must be from the dielectric such that energy is not coupled to the support structure. 
Second, and more important to this experiment, an understanding of the strength and profile of 
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the evanescent fields at the dielectric ends was created. The profile of the evanescent fields at 
the resonators ends determines the distance the conductor must be from the dielectric in order 
to couple energy to the conductor. 

Section 3.4 showed that placement of the conductor with respect to the dielectric had a 
dramatic effect on resonant frequency. This implied that the experimental structure must have 
had the ability to place different conductors at the same distance in a repeatable and accurate 
maimer. 

The last section yields insight into the goal of the NASA experiment, relating 
conductivity of the mesh conductor to the Q of the resonant system. Fig. 3.15 and 3.16 showed a 
weak relationship between the conductivity and the complex resonant frequency, suggesting only 
marginal usefulness of this measurement approach. In addition, the incidence angle on the 
conductor is not conducive to the desired near normal incident measurement. 

At this point, the work remaining is an iterative process of going from the model to 
experimental measurements, and then back to the model. In addition, accurate measurements of 
the dielectric constant, dimensions of the dielectric, and accurate measurement of air gap could 
be incorporated into the simulation to yield more accurate simulation results. To further refine 
both the simulation and to redesign the experiment, the dimension of the dielectric could be 
chosen to cause field propagation to be more planer. This would yield stronger evanescent fields 
and more coupling of energy to the conductor. This would also make the approximation of the 
TE boundary conditions at the dielectrics end's more accurate resulting in a better simulation. 
Calibration of the model should be done using conductors with known conductivity such as 


copper. 
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APPENDIX A 


DERIVATION OF TE COMPONENTS IN TERMS OF H z 


The following derivation provides the components and H p in terms of H f for the 
TE 0 , mode in the dielectric resonator of Fig. 1. The derivation begins with Maxwell's curl 
equations. 

VxE=-ju)iH (A.l) 

V x H = j <mE (A.2) 

If (A.l) and (A. 2) are expanded into cylindrical components, we obtain 


„ 1 dE, 

-j "» H P = l-§f-sr 

(A. 3a) 

iutE _ 1 dH * 
i E P ~ P &t> dt 

(A. 3b) 

„ dB » dE . 

(A. 3c) 

■ c dH o 

}“ E * = -d T~W 

(A.3d) 

- iwfiH, = J (jfypEJ - 

(A.3e) 


(A.3f) 


For the TE 01 mode, the following components are zero: 

E, = 0 E p = 0 H+ = 0 


Eqs. (A. 3) become 

dE d . QH 

j ^ H P~~ST ( A4& ) 0 = P ( A4b ) 
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0=0 


(A. 4c) 


c J H P M . 

jw< E * -~dT"di 


(A. id) 


M* H, = \fypE+) 


(A.4e) 



(A4f) 


Eq* (A. 4b) And (A.4f) imply that H fl and H t are not functions of <t>. Setting the <fi variation to a 
constant, Eqs. (A. 4) simplify to 


dE+ 

w H p * -dT 


(A.5a) 


jwpW, = 


( A .5b) 


and 


iUtE * ~~dT dp 


(A.5c) 


Taking the derivative of (A 5a) with respect to i and substituting into (A.5c) gives 



Each region must be considered separately in solving the form for (A.6). The resultant 
forms may then be combined to complete the problem solution. For regions, 1 and 2 of Fig. 1 we 

have 


^ (A.Ta) 

dr 3 *’ 


in region 3 we have 
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and In region 4 we have 


(A .7b) 



(A. 7c) 


We now develop the fields in each region for the TE 0l mode, beginning with region 1. 
Let ua define the magnetic field intensity as 


H] - A x sin *,(* - t 0 ) Jo(k\pP) > 


(A. 8a) 


for an arbitrary i 0 , where k lp is defined as 

- ** 

From (A.6) and (A.7a), we may write 

E i = ^ " O J i>(V> 


(A.8b) 


and from (A. 6a) 



-5T = MiV M 


- *„) jwm ■ 


In region 2 we define 

//* = A a sin *,(r - r 0 ) K 0 Q^) . 


(A.8c) 


(A. 9a) 
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where k^ p is defined aa 



Aa before we may write 


= jw/i 3 /t 3 i 3p jm k,{i - i 0 ) K„(j k 7p p) 


(A.9b) 


and 


H) = k t A^ p cot k t (: - t 0 ) Ki(i k 7p p) (A.9c) 

In region 3 we define (d = L) 


H 


3 

< 


= /t 3 e- T < 1 -^ J 0 (k 3pP ) , 


(A. 10a) 


where k ip is defined as 


k 3p = 


As before we may write 


= J'(^) 


(A. 10b) 


and 
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Wj= ’ 

[n the last region to be considered, region i t we define 


H] = A t sinh <(i + I) J 0 ( k ipp) > 


(A. 10c) 


(A.lla) 


where k ip is defined as 


k ip = + ? 

Again we write 


= jwfi^i^sinh <(i + l) Jq(* 4 ^) 


(A. lib) 


and 


- <^ 4 * 4 , C0,h <(' + 0 J o(* 4 ^) ' (A. tic) 

For simplicity of solution, a full modal expansion will not be done, but rather the fields 
in regions 6 and 6 will be neglected and continuity across a boundary will be applied as though 
the boundary were infinite. In other words, we set 

* 3 p ~ * 4 p ~ k \p- 

We have already imposed the equality condition on k t in regions 1 and 2. 

With these approximations, Eqs. (8a) through (11c) become the following: 
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Region 1: 


H\ - Jin *,(* ~ * 0 ) 


(A. 8a) 


E l = i"#*i ~dp 


JJ(V) 


H\ - r-*- ^ - k y A col *,(l - «„) 

' jwji, a* *ip 


//? = A 3 Jin t,(* - O K 0 (ji 3 ^) 


g» -Mhh, in *.(,-,„) Ki0*3^) 

* *3? 


4 ( x -^)Ki 0 V>- 

p hp 


h] = A 3 t ' d) J 0 {k lp p) , 


= j^) 


W 3 = _^J e -7(.-0 Ji(V ) 


Region 4: 


//} = A 4 finh <(x + 0 > 


(A. 11a) 
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SjJ-SJ^iiiiih <(* + *) y 0 ( k ^) 

9 *1 p 

and 

H 4 = cosh <(* + I) 

*U> 


where 

*lp = - k] , 

= \1 w3< jM 3 -”*] - 

> = \J*U “ . 

And 

< = >|^- W 4 « 4/‘ 4 • 


(A. lib) 


(A. 11c) 


(A. 12«) 


(A. 12b) 


(A. 12c) 


(A.12d) 
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APPENDIX B 

DERIVATION OF EIGENVALUE EQUATIONS OF SECTION 

2.2 

The boundary conditions for the resonator of section 2.2 are 

1) H](p=a) = H](p=a) 

2 ) B\(p=a) = El(p=a) 

3) //}(,=0) = H*(x=0) 

4) £}(*=0) = £}(i=0) 

5) H\(z=d) = Hfr=d) 

«) El(z=d) = El(z=d) . 

Applying each of these boundary conditions to the fields of App. A, we obtain the following 
equations: 

A, J 0 (V> = ^ K oGM ( B1 ‘) 
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^1 k Ip — A 2 p 

(B.lb) 

A x k f co* k t z 0 = A a < co»h < i 

(B.lc) 

- sin k t z 0 = j 4 4 finh (i 

(B.ld) 

A \ k , c °* *«( rf - 0 = - ^3 y 

(B.lc) 

A x sin k t (d - z c ) = A 3 

(B.lf) 


Since the eigenvalue* (in other word*, k t and 7 ) are de*ired in order to determine the 
resonant frequency, the amplitude coefficient* may be eliminated from the equation* to provide 
the simple eigen equations which must be solved. Dividing (B.lb) by (B.la), and recalling the 
relationship* of k Xp and k^ p to k t in (A. 12), we obtain an equation for k t in term* of the frequency 

w a* 




An additional equation is required if we are to determine either the resonant frequency 
or the k f . We may obtain such an equation by deleting and A ^ from the boundary equations 
also. Dividing (B.ld) by (B.lc) we obtain 


- tan k g i 0 = ~ tanh <J 


(B.3) 
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end dividing (B.lf) by (B.le) we obtain 


The last two equations have introduced the additional variable, of s, and y. The variable y is 
defined by (A. 12c). The unknown i 0 may be eliminated by combininc (B.3) and (B.4). 

To eliminate we first expand the tangent of (B.4) using a standard trigonometric 
identity to obtain 


tan 


tan (M) ~ (*,*,) _ *i 

k ,( d ~ V i + tan (k t d) tan (Vo) 7 ’ 


(B.5) 


or 


tan 


^ + tan (k g d) k 
1 - ^ tan (k,d) 


(B.«) 


Eq. (B.«) may be rearranged to give 

jr + £ coth <1 

t*n (k d) = ^ —7 

l-2jcoth<t 


or 


k t d = tan " ‘(^) + tan l (^ coth <l) 


(B T) 


Eqs. (B.2) and (B.7) form the two equations to be solved simultaneously for the 


eigenvalue of the system, k„ and thus the resonant frequency. 
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